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Abstract. A method is proposed of obtaining (2+l)-dimensional non-linear 
equations with non-analytic dispersion relations. Bilocal formalism is shown 
to make it possible to represent these equations in a form close to that for 
their counterparts in 1+1 dimensions. 
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1. Introduction 



Non-linear evolution equations with non-analytic (singular) dispersion rela- 
tions (SDR equations) form an important class of equations integrable by 
means of the inverse spectral transform. The Maxwell-Bloch equations [|TJ 
are the well-known representative of this class in 1+1 dimensions. A general 
construction of (l+l)-dimensional SDR equations solvable via the Zakharov- 
Shabat spectral problem was given by Leon [Q]. As was shown by Boiti et 
al in an interesting paper ||, the SDR equations in 2+1 dimensions pos- 
sess a number of peculiarities, the main one being the absence of an explicit 
expression for the evolution linear operator = dt — W which enters the 
Lax representation. Nevertheless, this circumstance does not prevent from a 
construction of soliton solutions by means of the Backlund transformations. 
In particular, it was proposed in || a (2+l)-dimensional generalisation of 
the Maxwell-Bloch equations which had a form of a rather complicated sys- 
tem of four equations. In our opinion, such a complexity was caused by the 
fact that the approach realised in was primarily based on the function 
W given unexplicitly. In this connection, it is seemed to be reasonable to 
propose another way of deriving the above class of equations without making 
the direct use of the function W. 

We will consider as a primary object a spectral transform R appearing 
in the framework of the 8- method [4-7]. Hence, the aim of the present pa- 
per is to obtain a hierarchy of (2+l)-dimensional non-linear equations with 
non-analytic dispersion relations compatible with the linear evolution of the 
spectral transform R. We will demonstrate that the formalism developed by 
Beals and Coifman || for holomorphic dispersion relations can be naturally 
adapted for equations of interest. Our consideration relies essentially on the 
bilocal approach initiated by Konopelchenko and Dubrovsky || and elabo- 
rated to a full extent by Fokas and Santini [|l(], [TT| . It is precisely the bilocal 
formalism that allows us to generate in a natural manner (2+l)-dimensional 
analogues of many structures which successfully work in 1+1 dimensions. 
We will show that the form of the SDR equations in 2+1 dimensions written 
in bilocal variables is very close to that for equations in 1+1 dimensions. 
In particular, the 'squared eigenf unction' structure typical for the (1+1)- 
dimensional situation takes place also in 2+1 dimensions. In the process of 
deriving a hierarchy of equations we shall not use, as distinct from [10|, an 
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extended integral representation for the function W (due to the lack of an 
explicit expression for it). In our view, the proposed way of obtaining the 
recursion operator follows more closely the lines of 1+1 dimensions. 

In the following, we shall restrict ourselves to the consideration of the 
hyperbolic spectral problem. A derivation of the relevant formulae in the 
case of the elliptic spectral problem does not cause principial difficulties. 

2. Lax representation and <9-problem 

As a starting point in a construction of nonlinear SDR equations, we consider 
a (9-problem on a complex plane C (3 = d/dk): 

B(f)(k) = J J dl Adl<f>(l)R(k,l), k,leC, (1) 
4>{k) = 1 + 0(1/*;), A;^oo, G SL(2, C). 

Here the matrix R (the spectral transform) is a distribution in C 2 and a time 
dependence is given by the following linear evolution equation: 

d t R{k, I) = R{k, l)Q(k) - Q(l)R(k, I). (2) 

In the above equation Q(k) is a matrix-valued function on C called a disper- 
sion relation. In a general case, Q(k) consists of a holomorphic (polynomial) 
part flp(k) and a non-analytic (singular) part fl s (k), i.e. dQ s ^ in a some 
subset of the plane C. 

Let us denote the integral in (1) as <p(k)R k F where F is an integral 
operator acting on the left in accordance with the equation (1). Hence, we 
write (1) as 

d<j>(k) = <p(k)R k F. (3) 

A solution of the 9-problem is given by a solution of the linear integral 
equation 

Mk) = 1 + -!- f [ d ] A f [ [ dm A dm <j)(m)R(l, m) 
2m J J I — k J J 

= l + ^-jj J^f W) R i F) = l+ <P{k)Rk FC k . 

The integral operator C k acting on the left transforms the argument 
function in front of it to I and integrates with the weight (27ri)~ 1 (Z — 
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the whole complex plane. The integral operators introduced in such a way 
allow us to write formally a solution of the <9-problem (3) as 

<j ) (k) = l-(l-R k FC k )-\ (5) 

The similar representation for solutions of the <9-problem was effectively used 
by Beals and Coifman || in the case of holomorphic functions Q p (k). 
Let us define a pairing for matrix-valued functions on C: 

(<f>, i/>) = — [ [ dk A dk 4>(k) $(k), 
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where tilde stands for transpose. With respect to this pairing we have 

(<f>C k , V) = -(0, VOk), (<PRkF, = (0, ^R k F), (6) 

where R(k, I) = R(l, k). Assume then a parametric dependence of R(k, I) 
on spatial variables (x, y) of the form 

d x R(k, I) = ila 3 R(k, I) — ikR(k, l)a 3 , 

d y R(k, I) = i(k - l)R(k, I). (7) 

Taking as a basis the representation (5), it is shown in Appendix 1 that 
the choice (7) is equivalent to setting the two-dimensional Zakharov-Shabat 
spectral problem 

7\0 = (d x + a 3 d y + Q)<f> - tk[a 3 , 0] = 0, (8) 

where a potential Q is defined as 

Q(x, y) = -i[a 3 , {(t>R k F)} (9) 

and we denote (/, 1) = (/). 

Now we turn to an evolution linear problem d t (j) = W<p + <fi Q. It follows 
from (5) and (2) that 

4> t = d t R k FC k (l - R k FC k y l = (<pR k FQC k -i> QR k FC k )(l - R k FC k y l 
= (0 J R fc FfiC fc -0fi)(l- J R fc FC fc )" 1 + 0fi, 
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which gives 

W<j>=(<j>SlR k FC h -<j>Sl)(l-R h FC h )- 1 . (10) 
It is shown in Appendix 2 that the formula (10) is reduced for Q p = to 

W(k)(f>(k) = -<f>(k)8n a (k)C k (l - R k FC k )-\ 

Multiplying this relation on the right by (1 — R k FC k ) and acting by the 
(^-operator, we obtain 

dW(f) + W((j)R k F) - W<j>Rk F = -<j>d O s , 

which gives the integral equation for the function W |?J: 



dW{k) = -<pdVL s <p-\k) + / / dlAdl [W{1) - W{k)} (f)(1) R(k, l)<p-\k). 

(11) 

Hence, the function W is known only to withih a solution of the integral 
equation (11). Nevertheless, Boiti et al have shown || that it is possible to 
derive SDR equations from the corresponding Lax representation with the 
operators T\ and T 2 = d t — W . 

It should be noted here that the equation (11) includes the inverse func- 
tion cf)^ 1 . However, there is no in 2+1 dimensions, contrary to 1+1, a simple 
equation (like (8)) for Hence, a problem arises to find a natural (2+1)- 
dimensional analogue of the inverse function in 1+1 dimensions. We shall 
show in the following section that such a function does exist and it permits 
to simplify the SDR equations in ||. 

3. Hierarchy and recursion operator 

Let us calculate the evolution of the potential Q which is given explicitly by 
(9): 

dtQ = -i[(T 3 , {d t (<pR k F))\. 
The right-hand side can be transformed as follows: 

dt (<f>R k F) = d t <f)R k F + cj) d t R k F = W(j>R k F + <\>R k FSl. 

Further calculation, due to (10), yields 

d t ((f>R k F) = <pR k FVtC k {l — R k FC k )~ l R k F — 00(1 — R k FC k )~ l R k F 
+4>R k FQ = <pR k FSl(l-C k R k F)- 1 -<pSlR k F(l-C k R k F)-\ 
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Hence, 

d t Q = -i[a 3 , (<f>R k Fn(l-C k R k F)-\ 1) - (<f)QR k F(l — C k R k F)~ 1 )]. 

Taking into account properties (6) of the pairing, we get 

d t Q = -i[az, ^RkF^l-^ + RkFCk)- 1 )-^^ 1 • (1 + R k FC k )~ l R k F)\. 

(12) 

Now we introduce a dual function <f>*(k) by means of the relation 

4>*(k) = l-(l + R k FC k )-\ (13) 
The (9-problem for the dual function has the form 

d(f)*(k) = - J J dlAdlR(l, k)(f)*(l), d<t>*(k) = -<j)*{k)R k F (14) 

and 4>*{k) satisfies the dual spectral problem 

d x <t>* + d y <f)*<j 3 - <j)*Q - ik[a 3 , 0*] = 0. (15) 

The derivation of (14) and (15) is given in Appendix 3. It is seen from (15) 
that just the dual function <fi* is a true (2+l)-dimensional generalisation of 
inverse functions in 1+1 dimensions. It should be stressed that the defini- 
tion (13) of the dual function arises naturally within the framework of the 
formalism based on the representation (5). 

Taking into account the above relations concerning the dual function, we 
write the evolution (12) in the form 

d t Q = -ifa, (<pR k Fn, 0*> - 0*R k F)] 

= -i[a 3 , (B<t> ^0*) + {<t>nB<t>*)\. 

Finally, dividing the dispersion relation into regular and singular parts, we 
obtain under condition Q s (k) — > for k — > oo: 

d t Q = -i[a 3 , («9(0fi p 0*)) - (48n s (f>*)]. (16) 

Assume further that 

flp(k) = a n k n a 3 , a n = const, n — 0, 1, . . . , 
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which yields dfl s (k) = g(k)a 3 . Now we introduce a bilocal object 

M 12 (x, y lt y 2 , k) = y x , k)a 3 (f)*(x, y 2 , k) = 0icr 3 02 

It is easy to see that the function M 12 satisfies the equation 

d x M 12 + a 3 d yi M l2 + d y2 M 12 a 3 - ik[a 3 , M l2 ] + Q 1 M 12 - M l2 Q 2 = 0, (17) 

where Qi = Q(x, y^), i = 1, 2. Hence, equation (16) takes the form 

Snd t Q 2 = -ta n 5 12 [a 3 , (8 {k n M 12 ))\ + iS 12 [a 3 , (g(k)M 12 )}, (18) 

where 5\ 2 = ${y\ ~ 1/2)- Following [[UJ, we introduce the notations 

P 12 M 12 = d x M 12 + a 3 d yi M 12 + d y2 M 12 a 3 , Q±M 12 = Q 1 M 12 ±M 12 Q 2 . (19) 

Let My 2 and Mf 2 are the diagonal and off-diagonal parts of the matrix M 12 , 
respectively. Then (17) and (19) yield 

P12 M^ + Q" M£ = 0, (20) 

P 12 M a l2 - 2 ika 3 M a l2 + Q l2 Mf 2 = 0. (21) 

We can write from (20) the diagonal part as Mf 2 = 03 — PviQxiMvi- Hence, 
the equation (21) is written in the form (A — k)M® 2 = (2i)~ 1 Qf 2 ■ 1, where 
the operator A is defined as 

A = ±a 3 (P l2 -Q u P 12 1 Q u ). 

Then Mf 2 = (2i)- 1 (A-k)- 1 Qf 2 -1 and after expansion (A-A;)" 1 = - £~ =1 fc^A" 1 " 1 
we can write the polynomial contribution to d t Q in (18) as 

00 

-ia n 8 12 [a 3 , (d(k n M 12 ))} = a n a 3 5 12 ]T (dk n ~ m ) A m ~ l Q+ 2 ■ 1 

771=1 

= - % -a n a 3 8 12 A n Q+ 2 ■ 1. 



Now we have all the needed to formulate a closed system of equations describ- 
ing the evolution of the potential Q under condition of the linear evolution 
of the spectral transform R: 

5i 2 d t Q2 = ~a n a 3 5 12 A n Qt 2 ■ l+i5 12 {a 3 , (g(k)M 12 )), 

(Pi2 M 12 - ik[a 3 , M 12 ) + M 12 ) g{k) = 0. (22) 

Here the operator A plays the role of a recursion operator (more precisely, 
A is connected with the true recursion operator by means of cr 3 [|10|)- If 



g(k) = 0, we get from (22) the well-known hierarchy including the Davey- 
Stewartson-1 equation derived by Santini and Fokas |]10| on the basis of an 
integral representation for W. 

For Q p = the system (22) takes the form 

5 12 d t Q 2 = i6 12 [a 3l (g(k)M 12 )}, 

(P 12 M 12 - ik[a 3} M 12 ] + Qr 2 M 12 ) g(k) = 0. (23) 

It is seen that the structure of the system (23) is similar to that for (1+1)- 
dimensional Maxwell-Bloch equations (g(k) ~ 5(Qk)5($lk — a)): 

d t Q = i[a 3 , (g(k)$(k))}, $ = 0a 3 0- 1 , 

d x §{a) - ia [a 3 , $(«)] + [Q, $(«)] = (24) 

and the system (23) is reduced to (24) in the (l+l)-dimensional limit. It 
should be stressed that, as distinct from [fj, the (2+l)-dimensional counter- 
part (23) of the Maxwell-Bloch equations demonstrates explicitly the pres- 
ence of the 'squared function' term. It can be shown that a function T 
introduced in is expressed, as a matter of fact, in terms of the above 
squared functions as T = id (kef) <$*) . Finally, for n = 2 equations (22) yield 
a (2+l)-dimensional generalisation of the equations derived in |T2| in the 
context of non-linear optics. 

4. Conclusion 

We proposed a procedure for obtaining (2+l)-dimensional non-linear equa- 
tions with non-analytic dispersion relations compatible with the linear evolu- 
tion of the spectral transform. An important step in deriving these equations 
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was to use the representation (5) of the <9-problem solution. In spite of the 
formality of this representation, it allows us to perform all the needed ma- 
nipulations. The introduction of the dual functions has made it possible to 
obtain the hierarchy of equations without explicit use of the second Lax oper- 
ator. The application of the bilocal formalism was crucial for bringing these 
equations to a form similar to that for their counterparts in 1+1 dimensions. 
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Appendix 1. Linear spectral problem 

We show here that the choice (7) of the dependence of R(k) on spatial vari- 
ables x and y leads to the Zakharov-Shabat problem on the plane. Differen- 
tiating (5) with respect to x we obtain d x <f> = (pd x R k FC k (l — Rk FCk)^ 1 . In 
virtue of the definitions of the integral operators F and C k and (7) we can 
perform the following calculation: 

(j)d x R k FC k = f ~JZ~j~ J J dm A drh(/)(m)d x R(l, m) 

= - — : J J — ^ J J dm A dm im (p(m)a 3 R(l, m) 

1 f f dl Adl , f f , _ . \ T-j/ 7 N 

— - — : J J — — —il J J dm A dm(p{m) R(l, m)a 3 

- h II *°> R ' F >-^. II d < A dU t 1 + rh) <* R ' F) - 

(A.l) 

Since we have from (4) <pR k FC k = — 1, then (A.l) and evident relation 
R k FC k (l - R k FC k )- 1 = {\-R k FC k )~ l - 1 yield 

d x <f) = -ik 0cr 3 -i((f)R k F}a 3 <j) + ik cr 3 (l - R k FCk)- 1 . (A.2) 

Similarly, 

d y (j) = ik(j) + i((j)R k F)<j> -ik(l- R k FC k y l . (A.3) 
Adding (A.2) and (A.3) yields 

d x <j> + a 3 d y <t> - ik [(7 3 , 4>] - i[a 3 , (<j)R k F)](f) = 0. (A.4) 

Hence, if we identify — i[a 3) (<f)R k F)] = Q(x, y), (A.4) gives the above spec- 
tral problem. 

Appendix 2. Linear evolution problem 

In order to derive the linear evolution problem d t (f) = W<f) + <f)Q, we calculate 
d t (f) from (2) and (5). Let us take for simplicity Q p = whereas 
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which gives dfl(k) = g(k)cr 3 . The calculation yields 

d t (t> = ct ) d t R k FC k {l-R k FC k y l 

= [<PR k FSlC k -<pSlR k FC k ](l - R k FC k y l 
= <f>R k FnC k (l-R k FCk)- 1 -<t>n{l-R k FC k y l + SI 

= w(k)(f)(k) + (f)(k)n(k), 

where 

W{k)4>{k) = ((PRkFnC.-tPtyil-R.FC,)- 1 . (B.l) 

Taking into account the definitions of the integral operators F and C k , we 
can rewrite (B.l) as 

W{k)<P{k)^-RkFC k ) = ^-[ 



2m J J I — k 



1 r f dl Adl f f , , _ , , , _ . 1 r f ds Ads . . , ^ 
= 2m J J T^k JJ dmAd ™ <K m ) R «> m ^JJ -J^M** - n. 

(B.2) 

The denominator in (B.2) can be represented as 



(l-k)(s-l) s-k\l-k l-sJ' 

Then we have 

W<j>(l-R k FC k ) 
1 r r dl Adl r r , . ,,,„,. , 1 f f ds Ads 



1 f f dl A dl f f . ,_ , , s x 1 f f ds A ds 

= 2mll T^k JJ dmA d ™«K™W> "») ^ y y t^f^ (s)(73 

1 f f ds Ads 1 r r dlAdl f f . ,_ . . . . , 

^ y y t^f^ (s) 2^ y y y y ^ a ^ ^ - ^ 

1 r r d ^ A d^ 

= 4>R k FC k n(k)-—J J ^^g(s) (<PR S FC s )a 3 -<f>n 



= -g(k)(j>(k)cr 3 C k , 
where, as in Appendix 1, we use 0i? s FC S = 0(s) — 1. Hence, 

W{k)<t>{k) = -g(k)</>(k)a 3 C k (l - R k FC k y\ 
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It should be noted that the calculation of dt4> for fl s = and Q 
on the basis of (5) leads to the well-known operator T2 = dt 
1 equation fl5L (□ 



p - a 2 k 2 a 3 
W for the 



Davey-Stewartson 
has the form 



where a potential A of the mean flow 



A = 2a 3 [(k{(f>R k F)) - (k(PR k F) -i{d y 0R k F)] c 



Appendix 3. Dual spectral problem 

Here we give a derivation of the equations (14) and (15). The definition 
(13) gives 4>* = 1 — (p*Rk FC k and taking into account the evident property 
df{k)C k = f(k) for any function f(k), this yields 

d4>* = -4>*R k F — — [ [dlA (U4>*(l)R(k, I) 



dlAdl4>*{l)R{l, k) 



dlAdlR(l, k)(j)*(l) 



(transpose) 

Hence, equation (14) follows. Now we find a spectral problem for the dual 
function 0*. Differentiating (13) with respect to x, we find d x <fi* = —cj)*d x R k FC k (l+ 
R k FCk)^ 1 . Taking into account that R(k, I) = R(l, k), we obtain from (7) 

d x R(k, I) = ik R(k, l)a 3 - il cr 3 R(k, I), d y R(k, I) = -i(k - l)R(k, I). 

Then following the calculation in Appendix 1, we obtain 

d x (f>* = i(f>*u 3 — i(4>*R k F, l)a 3 4>* -ika 3 (l + R k FC k )-\ 

and 

d x <P* = ik a 3( j)* - i0V 3 (l, 4>*R k F) -ik(l + R k FC k )7^,a 3 . 
Similarly, 

d y <P* = -ik<f>* + i<f>*(l, ^R k F)+ik(l + R k FC k ) { l iansp y 

Hence, 

d x 4>* + d y( j}*a 3 -ik[a 3 , (jf]+i(j>*[a 3 , (1, 4>*R k F)] = 0. (CI) 
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We need now a connection of (1, <fi*R k F) with Q. It can be found as follows: 

(<f>R k F, 1) = (1 • (1 - R k FC k y l R k F, 1) = (1 • (1 - R k FC k )-\ R k F) 

= (1, R k F(l + C k R k F)- 1 ) = (1, 1 • (1 + R k FC k )- l R k F) = (1, 4>*R k F). 
Hence, 

g = -i[a 3 , (<t>R k F, 1)] = -i[a 3 , (1, 0*^ fc F)], 
and we derive from (C.l) the equation (15). 
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